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Abstract We give evidence that the functional renormaliza- 
tion group (FRG), developed to study disordered systems, 
may provide a field theoretic description for the loop-erased 
random walk (LERW), allowing to compute its fractal di- 
mension in a systematic expansion in £ = 4 — d. Up to two 
loop, the FRG agrees with rigorous bounds, correctly repro- 
duces the leading logarithmic corrections at the upper crit- 
ical dimension d = 4, and compares well with numerical 
studies. We obtain the universal subleading logarithmic cor- 
rection in d = 4, which can be used as a further test of the 
conjecture. 



Keywords Loop-erased random walks 
Functional renormalization group 

PACS 05.40.Fb • 05.10.Cc • 64.60.av 



Sandpiles 



The loop-erased random walk (LERW) was introduced 
by Lawler JT| as an alternative to the self-avoiding walk 
(SAW), which is relevant in polymer physics. On the lattice, 
the LERW is defined as the trajectory of a random walk in 
which any loop is erased as soon as it is formed. As the SAW, 
the LERW has no self-intersections, but is more tractable 
mathematically. It has been proven that the LERW has a 
scaling limit in all dimensions [2 3 4 5 6 7]. The number 
of steps (or time) t it takes to reach the distance L scales as 
ti ~ U, where z is the fractal dimension of LERW. In d = 2 
this scaling limit is conformally invariant and described [8 
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9 1 by the stochastic Loewner evolution SLE2. Though the 
LERW and SAW belong to different universality classes, the 
LERW itself has received significant attention due to appli- 
cations in combinatorics, self-organized criticality (SOC), 
and, more recently, conformal field theory and SLE. The 
LERW can be viewed as a special case of the Laplacian ran- 
dom walk [ 10 1 and can be mapped [4] to the problem of uni- 
form spanning trees (UST): Chemical (i.e. shortest) paths on 
UST obey LERW statistics. It is proven that the upper crit- 
ical dimension is d uc = 4, the same as for the SAW, since 
for d > 4 the traces of two random walks do almost surely 
not intersect. Hence for d > 4 the fractal dimension of the 
LERW is that of a simple random walk, z = 2. It was proven 
a while ago [11] that for d < 4 it is bounded from above by 
the Flory estimate for the SAW exponent: 
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where R = {R{t) 2 ) l l 2 is the radius of gyration, and we have 
introduced £ = 4 — d. The mapping to UST and equivalently 
to the g-state Potts model at q — > was used in d = 2 to pre- 
dict 0] ZLERw(d = 2) = |, later proved in Ref. fl2). It is a 
particular case, for K = 2, of the fractal dimension df = 1 + f 
of the trace of SLE^- (a simple curve for K < 4). Another 
connection in d = 2 is to the 0(N) loop model with N = —2, 
both corresponding to a conformal field theory (CFT) with 
central charge c = —2. The leading logarithmic corrections 
at the upper critical dimension d = 4 where obtained by 
Lawler [5 ] who proved that: 

f(ln0 1/3 . 
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In d = 3 the value of zlerw is known only from numer- 
ics. The most precise estimate was obtained by Agrawal and 
Dhar fBI . 

zusw(d = 3) = 1-6183 ± 0.0004 , (3) 

improving on the previous numerical result 1141 z = 1 -623 ± 
0.011. 

Contrary to the SAW, which is described by the 0(N) 
model at = 0, there seems to be at present no field-theoretic 
approach to compute the LERW exponent in a dimensional 
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expansion around d = 4. This is surprising, especially when 
compared with the recent progress in CFT descriptions in 
d = 2. In this short paper we propose a field theoretic de- 
scription for the LERW, based on the Functional RG, a method 
developed to study disordered systems. We build on a con- 
nection proposed a while ago between the depinning tran- 
sition of periodic elastic systems, also called charge den- 
sity waves (CDW) in random media and sandpile models. 
The correspondence, on which we detail below, is indirect, 
through a chain of related models: From LERW to UST to 
sandpiles to CDW-depinning and finally to functional renor- 
malization group (FRG) field theory. Some of the connec- 
tions are not rigorous. At the end we show that the FRG 
passes all the tests of presently known results for LERW. In 
particular, it reproduces the correct leading logarithmic cor- 
rections in d = 4 given by Eq. (O, and makes a prediction for 
the subleading logarithmic correction which we hope will be 
tested in the near future. 

We now present briefly the intermediate models. The 
Bak-Tang-Wiesenfeld (BTW) sandpile model was proposed 
as a prototype for driven dissipative systems exhibiting SOC 
1151 . It is defined on the c?-dimensional hyper-cubic lattice 
with Lr sites. The configuration at time t is given by the in- 
teger number of grains h(x,t) at site x. The site x is unstable 
if h(x,t) > 2d in which case it relaxes according to the top- 
pling rule 

h(x,t + l) = h(x,t)-2d , 

h(y,t + l) =h(y,t) + l , (4) 
where y denotes all the 2d nearest neighbors of site x. The 
neighbor sites may then become unstable and the toppling 
continue. The order of topplings is irrelevant, thus the no- 
tion of "Abelian" sandpile, which allows to use e.g. paral- 
lel dynamics. The process continues until no unstable site 
remains, i.e. the avalanche ends. This is achieved through 
grains that leave the system at the boundary. To ensure a 
steady state one drives the system by adding a grain to a 
randomly chosen site x after each avalanche. Stable config- 
urations are either transient, or recurrent in which case they 
appear in the steady state with equal probability [ 16 1. 

Using the burning algorithm a one-to-one correspondence 
is set up between the recurrent configurations of sandpile 
models and spanning trees on the same lattice plus a sink 
xq (into which fall all grains leaving the system through the 
boundary) [17|. These spanning trees are connected sets of 
edges touching all lattice sites with no loops, and chosen 
with uniform probability (UST). Due to the Abelian prop- 
erty any avalanche initiated at an arbitrary site x, can be de- 
composed into successive toppling waves, such that during 
one wave each site topples at most once [18 |. Each wave is 
characterized by a two-root spanning tree on the extended 
lattice with connected jc; and xq. The subtree with the root at 
x, spans the sites toppled by the wave while the subtree with 
the root at xo connects the sites not affected by the wave. 
This spanning tree gives a unique (chemical) path between 
two toppled sites separated by the distance L. The number 
of edges along this path is nothing but the number of up- 
date steps between these two toppling events, i.e. the time t 



between them. Thus the time f and length L in the sandpile 
dynamics are related by the fractal dimension of the chem- 
ical path along a spanning tree [19], which is nothing but a 
LERW. 

Narayan and Middleton proposed that the charge density 
wave (CDW) near the depinning transition can be viewed 
to some extent as a BTW sandpile model H20II211 . The con- 
figuration of a fi-dimensional elastic object, such as a CDW, 
pinned by random impurities is parameterized by a displace- 
ment field 6 Mr. Here we assume that the distortions 
of the object caused by disorder are continuous. For CDWs 
u xt is the local phase divided by 2n. In the continuous limit, 
the dynamics is described by the over-damped equation of 
motion 

J]d,u xt = V 2 u xt +F(x,u xt )+f, (5) 

where V\ is the friction, / the driving force, e.g. the exter- 
nal electric field for CDW, and F the impurity pinning force. 
The latter is taken to be Gaussian with zero mean and corre- 
lator 

F{x,u)F{x' 1 u l ) = A{u-u')8 d {x-x'), (6) 

where A (u) is an even periodic function (chosen here with 
period 1). For CDWs F(x, u) is often written as 

F{x,u) = a(x)Y[u-y(x)] (7) 

with \Y(u) \ < 1 being aperiodic function with period 1. The 
strength a (x) of impurities and their selected phase y(x) £ 
[0, 1) are taken to be random and uncorrelated in space. Both 
forms (O and © of disorder correlations are equivalent at 
least in the vicinity of the so-called depinning transition. 
The system driven by a homogenous force / undergoes a 
depinning transition at a critical force f c : below the transi- 
tion / < f c the system is pinned while above it, it exhibits a 
never-ending motion. This dynamic transition is reminiscent 
of an ordinary second-order phase transition with the veloc- 
ity playing the role of the order parameter: v ~ (/— f c )P 
and enjoys a similar universality, i.e. independence of the 
scaling behavior on the microscopic details. The correla- 
tion length 4 defined through the velocity-velocity correla- 
tion function diverges at the transition as £, ~ (/ — f c )~ v \ 
the length and time scales of fluctuations are related by the 
dynamic critical exponent z: t ~ x z . Discretizing space and 
time, one can rewrite Eq. (O near the depinning transition 
in the form of the automaton model <(4j with h(x,t) playing 
the role of a coarse-grained curvature of the elastic object. 
For a pinning potential with narrow and steep wells the con- 
figuration of the CDW can be encoded by integer variables 
m(x,t) — u(x,t) — y(x). Upon time and space discretization 
m(x,t) evolves according to the following dynamic rule II20II 

m{x,t + \) =m(x 7 t) + e{f+V 2 [Y(x)+m(x,t)] + a(x)},(8) 

where © (x) is the Heaviside step function and V 2 is the dis- 
crete Laplace operator. In terms of a coarse-grained curva- 
ture defined as 

*(*) = £[m(y) -m(x)} +int[/ + V 2 y(x) + a{x)\, (9) 

y 
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Table 1 Pade approximants for z in d = 3. 
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the dynamics rule (HJl can be expressed in the form of the 
toppling rule <(4j with an inessential difference that the sta- 
bility threshold now is h(x) — not h(x) = 2d. Increasing 
/ toward the critical force is equivalent to increasing h by 
1 at a random site. Thus both models, the original BTW 
model and the discretized version of Eq. (01, are driven by 
adding grains, h(x,t + 1) = h(x,t) + 1. In the CDW model 
grains are added with a cycle restriction: a second grain can 
be added to a particular site only if all other sites have re- 
ceived a grain. The ordering of sites in the cycle is provided 
by the value of the quenched random term V 2 y(x) + a(x) in 
Eq. (|9). While this difference between the pinned at thresh- 
old CDW and BTW may seem inconsequential for large sys- 
tems, it illustrates that the mapping is presently not rigorous; 
moreover an ad-hoc discretization was used. Nevertheless it 
is supported by numerics H201 . It strongly suggests that the 
dynamic exponent z describing the depinning transition of a 
fi?-dimensional CDW coincides with the fractal dimension of 
LERW in d dimensions. 

The field theory which describes the depinning transition 
is based on the Functional RG H221 . Applying the Martin- 
Siggia-Rose formalism to Eq. (0! and averaging over disor- 
der with the help of Eq. © one arrives at the action 



S= ii2 xt (rid t -V 2 )u x t- / iu xt f xt 

Jxt Jxt 

iu xt iu xt iA(u xt - u xt i), 
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where the response field u xt has been introduced. We ar- 
gue that the field theory defined by ( TTOb also describes the 
LERW. The dynamic exponent of model d 1 01 ) at criticality 
is nothing but the fractal dimension of LERW. We also ex- 
pect that more properties of LERWs are encoded in the fixed 
point A*(u) since it controls the large-scale behavior of the 
model. At this stage, however, we only extract the relation 
between time f and space L, more work is needed to eluci- 
date relations to other observables. Power counting shows 
that the upper critical dimension of the model is d uc = 4. 
The peculiarity of the problem is that all derivatives of A (u) 
at the origin u = become relevant operators below d ac = 4. 
Thus, instead of a finite number of coupling constants we 
have to deal with a whole coupling function A (u). To extract 
the scaling behavior one has to study the flow of the renor- 
malized function A(u) under changing the infrared cutoff 
towards infinity. Recent progress has shown that full consis- 
tency requires a 2-loop study. The corresponding flow equa- 
tions to 2-loop order read 1231 

d 1 A{u) = eA{u)- l -[(A{u)~A{0)) 2 ]" 



-~[(^(«) 



-A(0))A'(u) 2 }"+ 1 -[A'(0 + )} 2 A"( u ), 



(11) 



d,lnr] = -A"(0)+A"{0) 2 + A'"{0 + )A'(0- 1 



ln2 



(12) 

where / = lni, and L the infrared cutoff, e.g. the size of the 
system. Below d = 4 the flow equation dill) has a fixed point 
(FP) solution A*(u) with a cusp at the origin: A*'(Q + ) ^ 0. 
Taking into account that f ~ T]/L 2 and that in the vicinity of 
the FP T]i scales with L according to Eq. d 1 2b . one finds that 
f ~ U, and to 2-loop order the exponent z is given by 

,2 



e__e 1 
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(13) 



Comparison with the exact bound OJ is encouraging since 
the 2-loop correction has the correct sign. To estimate val- 
ues of z in d = 2,3 we compute different Pade approximants 
[n/m] to (O. In d = 2 we obtain z(d = 2) = 1.23 ±0.22, 
consistent with the presumed exact value ^ but of poor ac- 
curacy due to the large expansion parameter e = 2. Ta- 
ble Q] contains the approximants [n/m] in d = 3. To improve 
the accuracy we construct the higher order Pade approxi- 
mants with n + m = 3 imposing that z(d = 2) = | for e = 2. 
The average and root-mean-square of all approximants with 
n + in = 3 yields 

z= 1.614 ± 0.011 . (14) 

The same procedure based on one loop only produces z = 
1.638 ± 0.012. The value (O is our best 2-loop FRG pre- 
diction for d = 3 and is in fairly good agreement with the 
numerical result OJ. 

The field theory is most predictive at the upper criti- 
cal dimension, where it yields exact results. For FRG some 
were obtained previously (see e.g. Refs. [24 25 26]). Here 
we compute the logarithmic corrections to the dynamics which 
yield a prediction for the LERW in d = 4. For £ = one 
shows from Jilt that the (periodic) disorder correlator ap- 
proaches the FP solution A * (u) = as follows (for < u < 
1): 



Ai(u) = 
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Substituting Eq. ( fT3T l into Eq. ( fl2l we obtain 
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Renormalizing the relation / 
at 

i , x 1 n r 21nlnL 



TJ/L up to scale L we arrive 
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which can be rewritten as 
lnlnf 
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(17) 



(18) 



1 We used the average and root mean squared of the set zn/o] = 4/3, 
Z[o/i] = 3/2, z[2/o] = 8/9, Z[ /2] = 6/5- We excluded approximant [1/1] 



which accidentally vanishes. 
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The scale I? can be taken as R 2 , the radius of gyration. We 
note that the leading order of Eq. ( 118b coincides with the 
result (|2]l of Lawler. Here we obtain the universal subleading 
correction. 

The prediction d 1 8b could be tested numerically, very 
much as for the corresponding prediction 11271 for the SAW 
in d = 4, checked in [28]; there subleading corrections are 
necessary to properly fit the numerical data at any feasible 
chain length. To this purpose, it is useful to note that, as for 
the SAW, there is only a single non-universal constant in the 
correction term in the parenthesis in ( 118b , c/lnf. It can e.g. 
be put to zero by a proper choice of to, setting f — > f/?o (in 
which case further corrections are 1 / In 2 1 and determined by 
3- and higher-loop terms, not considered here). For compar- 
ison, one similarly finds, from the 2-loop /3 -function 11291 
of the 0(N) model in d = 4 that L 2 ~ f(lnf) aiV (l - b N 1 -^-) 
with a N = [N + 2)/(N + S) and b N = (N + 2) (68 + SN - 
N 2 )/(N+S) 3 from which the Duplantier values for the SAW 
1271 , a Q = 1/4 and b = 17/64 are retrieved at N = 0. Note 
that no value of N can account for ( 1181 . thus a representa- 
tion of LERW via the 0(N) field theory, if feasible at all, 
would at least require a more complicated operator corre- 
spondance. 

Note finally that the exponent t for the avalanche-size 
distribution^ was recently computed to one loop within FRG 
13011311 . In the sandpile literature no controlled calculation 
exists for the corresponding exponent, usually called T s , but 
the formula X s (d) = 2 — 2/d leading to t s = 4/3 in d = 3 
has been conjectured [13] from scaling arguments. The FRG 
was found to agree to 0(e) with this formula, assuming T = 
T s . Evaluation of 2-loop corrections is in progress as a further 
test of the conjecture and of the relations between sandpile 
models and depinning. 

To conclude, it may appear surprising that LERW be re- 
lated to a field theory based on functional RG, where one 
would expect a "simpler", more conventional, field theory 
based on a single relevant coupling constant, as for the SAW. 
One clue may be that for periodic systems (CDW) the FRG 
posesses a stable submanifold with only two coupling con- 
stants, A(u) = a + bu(l — u) (for < u < 1), which con- 
tains the leading critical behavior. From (II lb one obtains 
1231131 1 the spectrum of convergence to the fixed point as 
l/(lnL) 1+a " in d = 4 and L^ m " with (O n = a n e + /3„£ 2 + 
0(e 2 ) fore/ < 4; a n = ±(3+n)(l+2«), ]8 n = -\{n+2)(2n+ 
1)(2« + 3). The values n = 1,2,... correspond to the con- 
vergence to the submanifold, and n = to the convergence 
inside it. Convergence to the submanifold is fast, with lead- 
ing eigenvalue 1 + (X\ = 5. Fast convergence was also found 
in numerics [32] where the FRG function A(u) in (fTTT l was 
directly measured in d = 1 , 2, 3. The conjecture raises many 
other interesting issues to be explored, such as the possibil- 
ity to predict other LERW observables and corrections to 
scaling, the role of other universality classes for CDW and 



2 The probability to have an avalanche of size s is proportional to 
s~ z , with a cutoff for large and small avalanches. 

3 Apart from the trivial uniform mode of eigenvalue e. 



interface depinning, the connections to fermionic field the- 
ory (known for UST, see e.g. [33], and conjectured for the 
FRG i26ll "). Work is in progress in these directions. 
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